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The angular momentum (partial wave) reduction of the Lippmann-Schwinger 
equation describing the interaction of two spin ^ particles is extended to the case in 
which the spin singlet and triplet states are coupled. A straight forward method for 
obtaining the angular momentum decomposition of the general potential is indicated. 
The derived formalism is needed to describe the interaction between two nonidentical 
spin-^ particles or between two nucleons when isospin symmetry is violated. The 
resulting modification of the Stapp phase-shift analysis is given. 
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I. INTRODUCTION 



Recent experimental advances are permitting detailed measurements and analyses of the 
scattering of polarized nucleons from polarized nuclei [|TJ. This, in turn, is leading to a rig- 
orous test of the theory describing spin-spin scattering 0. The general form of the elastic 
scattering T matrix for two spin ^ particles can be quite complicated. It simplifies consid- 
erably if we assume rotation invariance, parity conservation, and time reversal invariance, 
in which case 0-0: 

2T(k>, k) = a(k>, k) + b(k>, k) + (a- &K<r* + (c(A?, k) + d(k', k)) B> m & m 

+ (c - d)a\a\ + e(k\ k)(3> + 4) + /(£, k){<P n - <?*) (1) 

Here k and k' are the initial and final momenta in the center-of-mass system, <j p and a 1 are 
the Pauli spin operators for the projectile and target respectively, and the subscript on each 
a indicates a dot products with one of the three independent unit vectors: 

U v h 1 h 1 h ^ h -J- h 1 

n = — — , m = — — , / = — — (2) 

\kxk'\ \k'-k\ \k + k'\ 

For example, <t^ = a l -m. For experimental (on-energy-shell) scattering, the a—f coefficients 
in (JI]) are complex functions of the energy and scattering angle. Once these coefficients are 
known, all 36 possible experimental spin observables can be calculated from them []3],|| . 

Just as (d) is the most general scattering amplitude expressed as an operator in the direct 
product spin spaces of the projectile and target, so the projectile-target potential must have 
the same form: 

2V(k', k) = V a+b (k', k) + V a _ b (k', k)a p n a l n + V c+d (k', k)ala l m 

+V c _ d (k', k)afaj + V e (k>, k)(<% + a l n ) + V f (k', k){a p n - a l n ) (3) 

where the subscripts on V indicate the correspondence with ([[]). If the potential has the 
structure (|), then its use in a wave equation will generate a scattering amplitude of the 
form (0). In traditional language, the V a+ b term is a central potential, the V a _b, V c +d, and 
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V c -d terms are tensor forces (dipole-dipole interactions for two electrons), the V e term is 
the usual spin-orbit potential ||, and the Vf term is the unusual spin-orbit potential which 
couples spin singlet and triplet states. 

Prior angular momentum (partial wave) analyses, such as those found in Goldberger and 
Watson g, Stapp [g] and Goddard ||, have examined the scattering of two spin— | particles 
for the nucleon-nucleon case where the / amplitude vanishes due to the isospin symmetry. 
In our work on proton scattering from spin— | nuclei, we have had to develop a Stapp-like 
partial wave analysis JlO| for cases in which the / term does not vanish (extensions of the 



analysis of observables including an / amplitude had already been given by La France and 
Winternitz While the work of Gersten [[/J provides a partial wave analysis, it is in 

the helicity — not angular momentum — representation. While we suspect that others have 
deduced the general angular momentum analysis in the course of their research, we have 
been unable to find those analyses in the literature and so wish to present ours. Explicitly, 
we shall derive the partial-wave decomposition of the Lippmann-Schwinger equation when 
there is mixing between the spin singlet and triplet states (as well as within the triplet), and 
indicate how we have applied the method in our calculations. 

II. ANGULAR MOMENTUM ANALYSIS WITH SINGLET-TRIPLET MIXING 

Many-body effects and relativity leads to a potential V incorporating complicated non- 
localities. For this reason, rather than solve a Schrodinger equation we solve the equivalent, 
momentum space Lippmann-Schwinger equation: 

T(k', k) = V(k', k)+J E+ * P V(k',p) T(p, k) (4) 

where E{p) = E p (p) + E t {p) is the projectile plus target energy and the + superscript 
indicates a positive ie has been added to the on-shell energy E = E p (k Q ) + E t (k ). Since 
the solution of a three-dimensional integral equation is a bit too formidable, we reduce it to 
coupled one-dimensional equations by expanding T and V in partial waves fill: 
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v(k',k)=* £ ^v^' s \k',k)y^(k')yjr(k) (5) 

jrujll' ss' 

T{k>,k) = 2 - J2 ^- l) nt\k',k)y^(k>)yJrik) (6) 

Here / and s are the total orbital angular momentum and total spin for the target plus 
projectile, and j is the total angular momentum: 

j = l + s (7) 
i?= l( a l + a p) ; s = (s),l(t) (8) 

At times we denote the s = 0, singlet state with "s" and the s = 1, triplet state with "t". 
The y's in are spin-angle functions and are given in our conventions [|TT|] as: 

yt(k) = £ (lmism s \jm)Yr i (6,<P) \sm a ) (9) 

m s mi 



(2/ + l)(i-m)! qs 
\ 4vr(/ + m)! ' v ' K ' 



^(x) = (i-x 2 r /2 ^^ (ii) 

To express the Lippmann-Schwinger equation (f|) in the partial wave basis, we substitute 
(H) and (Q), the partial-wave expansions of T and V, into (f|) and employ the orthogonality 
of the spin-angle functions, 

/ yjr(k) yi%'(k) dk = 5 w 8 ss ,5 mm , (12) 

There results: 

nr\k\k) = V^\k',k) + lY, jyd P V ^^ p ^^ k) (13) 



LS 



In equation (|13D, the sum over L S is over the orbital angular momenta and spin states which 
are allowed to couple for a fixed value of the total angular momentum j. Accordingly, the 
integral equations to solve in the partial wave basis are (we leave off the (k', k) dependence 
of the leftmost T's and V's): 
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yj( ss ) 
v jj 
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jj 
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jj 
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jj 
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3 3 
3 (si) 
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(15) 



+ 



p 2 dp 
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+ 



p 2 dp 



E+ - E(p) 
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3+1 3+ 



i(A/,p) V]W_ x {k',p) 
V^f j+1 (k',p) V^)_ x {k' lP ) 



T 



j+ij+i 



(p,k) 



>j(tt) 



j-i j+i(P> k) 
(17) 



Equations (|I^-(|T5|) describe spin singlet-triplet coupling arising from the Vf term in the po- 
tential @ which in turn will produce an / amplitude in T ([[]). Equations (pp-fllTD describe 
coupling within the spin triplet state arising from the tensor force terms V a -b, V c+ d, V c -d in 
the potential (H), it mixes orbital angular momentum states |T2|]. Since the total angular 
momentum j is a conserved quantity, all coupled states have the same j superscript. 



III. EVALUATION OF V L ) L ' 

We use a two-step procedure to determine the potential matrix elements vff s \k', k) 
needed in (|I~4|)-(|TT|). First we evaluate the potential (|j) in the spin basis \s, m s ) and then we 
invert the angular momentum decomposition of these spin-basis potentials. The spin matrix 
elements required for the V a — V e terms of ([![) are given in Tables 7. 1-7 A of Goldberger and 
Watson U . As expected, the V a — V e terms have vanishing matrix elements between singlet 
and triplet states. The new term Vf has nonvanishing matrix elements only between singlet 
and triplet states, explicitly 
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(0,0\(< -#)IM> =(0,0|«-a*)|l,-l) =-<v^ (18) 

There is no 9 or dependence in (^Sj) because we adopt the "Madison Convention" shown 
in Figure 1. This convention takes the z-axis as the beam direction k (fa — 9% — 0), and 
places the scattered momentum k' in the xz plane (Of = 9,<pf = 0). Accordingly: 

v ss (k',k) = (o,o|v|o,o) 

= V a+b (k>, k) - V a . b (k', k) - V c+d (k', k) - V c . d (k', k) (19) 
V sl (k\k) = (0,0\V\1 A) = -V ls (k',k) = ^.,(9^) = ^V f (k',k) (20) 

v 00 (k', fc) = (i,o|y|i,o) 

= V a+b (k', k) + V a „ b (k>, k) + (V c+d (k>, k) - V c „ d (k>, k)) cos 9 (21) 
V u (k>, k) = VL X _! = V a+b (k', k) + V c+d (k', k) sin 2 °- + V c . d (k', k) cos 2 °- (22) 



V 10 (k',k) = ~V^ 



10 



■^V e (k>, k) - 4= k) sin9 + -^V c . d (k', k) sintf (23) 



V Q1 (k',k) = -V . 



i 



= ^=V e (k', k) - -^V c+d (k', k) sinfl + -^Vc-d(k', k) sin 5 (24) 
fc) = K.u = -K- 6 (^', fe) + fc) cos 2 ^ + V c _ d (&', fc) sin 2 °- (25) 

Here we have the matrix elements of V in the spin basis. We next expand these ma- 
trix elements in angular momentum states in order to determine the partial-wave matrix 
elements V^^ s \ We take the expansion of the potential in spin-angle functions @ and 
evaluate the matrix element between spin states: 

O 

/ /// WIT/V7/ 7M // "\ / HI HI I ' '\ 'I— I' 111 II l\- \ 

(s m s \V (k, k)\s m s ) = — \ s m s \sm s )i (I m z s m s \jm) 

^ j s' sl'lm s m g imim t i 

x (jm\lm lS rn s ) Yp' (9 f , 4> f ) Vtf'(V, k) Y™* (9 h fa) (sm s \s", m" s ) (26) 

where (k, 9i, fa) and (k 1 , 9f, <fif) are the spherical coordinates of the initial and final momenta. 
The Clebsch-Gordon coefficients vanish unless I, V = j ± l,j and rrij = mi + m s = m[ + m' s . 
Parity conservation requires |/ — l'\ = 0,2. In the "Madison Convention" (Figure 1) the 
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projectile has no angular momentum in its propagation direction and so mi = 0, in which 

case 



lf><^)=^(0.0) = y^ (27) 

As a sample, we concentrate on the new term V s \ ( p0|) which couples the 1 00) singlet 
state to the |11) triplet state. Because j is a constant and s' = in the final state, the total 
angular momentum j must equal The parity constraint then requires that 1 = 1'. Because 
m s and mi are 1 and 0, respectively, we deduce that rrij = 1, and my = 1. For this V s \ term, 



the sum in (p6|) reduces to a simple sum in the final orbital angular momentum 



E .••-£••■ (28) 

js' sll'm B m B imimy I 

Similarly, when the spin-dependent singlet- singlet potential V ss (k',k) and those within the 
triplet state V m > m (k', k) are evaluated, we obtain the desired partial-wave expansion of V 
(and of T with the interchange V — > T): 

V sl (k',k) = Z^^P/(x = cobQm) 2 l±L v i f t \k',k) (29) 
47r i=i # + 1) 

V ss (k\ k) = -^J2 Pl(*)W + 1) *0 (30) 

= ^E^) {(l + 2)V l l l " m (k\ k) - ^(l + l)(l + 2)V l l + + T\k',k) (31) 
+ (2/ + WT\k\ k) + (l- W l - m (k>, k) - V^Tik', k)} 

v 00 (k>,k) = ^ 2 Y,Pi(x){(i + i)v^ litt \k' 1 k) + iv l l l - m (k',k) 

Z7V 1=0 

+^(i + m + 2)v l %f t \k\k) + y/(i-i)ivtrT\k>,k)} (32) 

v 10 (k>,k) = ^Y.Pl^){-vl l - 1{tt \k\k) + v^ tt \k\k) 

47r 1=1 

+^v;+r\k', k) -^v;;_T\k\ k)\ (33) 

k) = ^ E P/(x) {-^T^V, fc) + fc) (34) 



1 

47T' 



1=2 



1 



X) 



l + l 



-V, 



i+i{tt) ( ,, 



II 



(k',k) 



1 



-V!tT\k\k) 



_ (Z + l)(Z + 2) 
V, Jfc) + Iv^T^CV, k) - J=fu-T\k>, k) 



(35) 



Note that the sum is actually over the orbital angular momentum V of the final state, 
but for notational simplicity we have changed the label V to I. Furthermore, note that 
the organization in (]29"|)-(|3~5|) combines matrix elements which multiply the same Legendre 
polynomial even though the matrix elements may correspond to different j values. The 
Lippmann-Schwinger equations (|I4|)-(|iT|) of course only couple states with the same j. 

We invert equations (|29])-(|3~5|) for the partial wave potentials V^ 8 (k', k) by projections 
based on the orthogonality of the associated Legendre polynomials. Specifically, we multiply 
the equation for each V m i m by Pj m ~ m ', and evaluate numerically the integral 



Im'm(k , k) 



-1 



dx V m > m (k', k) Pj m m| (cos Qy 



(36) 



For (|^) and (|30"D the inversion is simple because only one Vyrf ^ is involved: 

vT\k\ k) = v;t\k\ k) = -^kLi sl (k>, k) 



V l f s \k',k) = n 2 I ss (k',k) 



(37) 
(38) 



The equations fl3"T|)-(|3"5]) contain VJ?j 's intermixed for differing j and I values, and so the 
projection results in five coupled equations in five unknowns: 



I = BV 



(39) 



hi(k',k) 










I 00 (k',k) 








v!! tt \k',k) 


ho(k',k) 




TDm'm 




v!r m (k',k) 
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(40) 
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where [-B^ m ] is the matrix of coefficients multiplying the V^'s in (|3l|)-(j35|). In our application 
13| , the matrix equation ( f40l) is solved by numerical inversion 



V = B~ L I (41) 



and then the whole procedure is checked by recombining the potential according to (|29|) -(|35"1). 

Once the Vjf/ S ^'s are known, the coupled integral equations (fj4])-(|T7|) are solved numer- 
ically for the lf,\ B s ^'s using a modification of the Haftel-Tabakin technique |12| . Once the 
Ti,\ s s ^'s are known, the matrix elements in the spin basis {s'm s >\T\sm s ) are computed via 



the generalization of equations (|29|)-(|35D in which the V's are replaced by T's. Finally, once 



the (s'm s i\T\sm s ) are known, the a — f amplitudes of ([D) are obtained 0: 

a(k>, k) = ± (T n (k', k) + Taofi, k) - T^k', k)) (42) 

b(k', k) = l - (T u (k', k) + T ss (k>, k) + T^k', k)) (43) 

c(k', k) = l - (T u {k>, k) - T ss (k>, k) + T x _ x {k>, k)) (44) 

d(9,k) = \ (T 00 (k',k) + T 1 _ 1 (fc',fc) -T n (k',k)) /2cos^ fc (45) 

k) + T 01 (#, k)) /V2sin^ fc (46) 



e(k>, k) = (r 10 (#, fc) - T 01 (*', fc)) (47) 
f(k\k)=tV2T sl (k',k) (48) 

All spin observables are then calculated from a — f using the relations found in La France 
and Winternitz IHI. 



IV. RELATION TO PHASE SHIFTS 

The on-energy-shell T matrix elements in the partial wave basis T^j s s \ko, ko) can be 
related to phase shifts, a convenient phenomenological parameterization of the scattering 
data. This is particularly useful when the present formalism is applied to the two nucleon 
problem because there are tables of NN phase shifts. The relations to the bar phases are 

BIB 



-2ipT^ ss \k ,k ) = cos27 Z e 2 ^-l (49) 

-2ipTlf\k ,k ) = cos27 /e 2 ^-l (50) 

-2ipT^f j±1 (k ,k ) = cos2^e 2 ^' ±1 ^-l (51) 

-2ipT^ m (k ,k ) = -isin2^e 2i ^- 1 ^+ 1 ^ (52) 

-2zpT_/'f ) (fc ,fc ) = -ian^e^ 4 ^ (53) 

p = (54) 

^ °E p {k ) + E t {k ) 1 ; 

The parameter 7/ is the mixing angle between the 1 0, 0) singlet and |1, 1) triplet state, and 
the parameter e; is the mixing angle between the I and I + 2 triplet states. In Table I we 
give the connection to the a notation of Stapp || as well as that to the N spin notation used 
in our computer code Lpotp2. 

V. SUMMARY 

We have extended the partial wave analysis of the Schrodinger equation describing the 
interaction of two spin | particles to the case where a potential couples the spin singlet and 
triplet states as well as coupling within the triplet state. In particular, we have concentrated 
on the scattering configuration described by the Lippmann-Schwinger equation although the 
same formalism can be used for bound states. While a previous formalism was appropriate 
to the nucleon-nucleon problem, extensions are necessary when the two Fermions are not 
identical. We are now applying the new formalism to describe polarized proton scattering 
from polarized 3 He and 13 C nuclei |14|]. Our formalism should also be useful in describing 
the interaction between identical hadrons when isospin symmetry is violated. Furthermore, 
the formalism may find some applicability in the electron-atom interaction in atoms such as 
calcium where there is strong singlet-triplet mixing. 
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FIGURES 

FIG. 1. The coordinate system used to describe the scattering of momentum k into k' in the 
Madison convention. The incident momentum is along the z axis and the final momentum is in 
the xz plane. Note that 6' is the same as the O^k i n equations (p29j)-(|35|). 
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TABLES 
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